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Information about the neutron-star equation of state is encoded in the waveform of a black hole- 
neutron star system through tidal interactions and the possible tidal disruption of the neutron star. 
During the inspiral this information depends on the tidal deformability A of the neutron star, and we 
find that A is the best measured parameter during the merger and ringdown as well. We performed 
134 simulations where we systematically varied the equation of state as well as the mass ratio, 
neutron star mass, and aligned spin of the black hole. Using these simulations we have developed an 
analytic representation of the full inspiral-merger-ringdown waveform calibrated to these numerical 
waveforms, and we use this analytic waveform to estimate the accuracy to which A can be measured 
with gravitational-wave detectors. We find that although the inspiral tidal signal is small, coherently 
combining this signal with the merger-ringdown matter effect improves the measurability of A by 
a factor of ~ 3 over using just the merger-ringdown matter effect alone. However, incorporating 
correlations between all the waveform parameters then decreases the measurability of A by a factor 
of ~ 3. The uncertainty in A increases with the mass ratio, but decreases as the black hole spin 
increases. Overall, a single Advanced LIGO detector can measure A for mass ratios Q = 2-5, black 
hole spins Jbh/Mb H = —0.5-0.75, neutron star masses Mns = 1.2Mg— 1.45Mq, and an optimally 
oriented distance of 100 Mpc to a 1-a uncertainty of ~ 10%-100%. For the proposed Einstein 
Telescope, the uncertainty in A is an order of magnitude smaller. 

PACS numbers: 97.60. Jd, 26.60.Kp, 95.85.Sz 



I. INTRODUCTION 

By the end of the decade a network of second genera- 
tion gravitational-wave detectors, including the two Ad- 
vanced LIGO (aLIGO) detectors [TJ, Advanced Virgo [2], 
KAGRA [3] (formerly LCGT), and possibly LIGO- 
India 4], will likely be making routine detections. Fu- 
ture ground based detectors such as the third generation 
Einstein Telescope (ET) [5] , with an order of magnitude 
higher sensitivity, are also in the planning stages, and 
may be operational in the next decade. A primary goal of 
these detectors is extracting from the gravitational wave- 
form information about the sources. Of particular inter- 
est are compact binaries whose waveform encodes the 
sky location, orientation, distance, masses, spins, and for 
compact binaries containing neutron stars (NS), informa- 
tion about the neutron-star equation of state (EOS). 

The study of EOS effects during binary inspiral has 
focused mainly on binary neutron star (BNS) systems. 
Work by showed that EOS information could be im- 
printed in the gravitational waveform through tidal inter- 
actions. In the adiabatic approximation, the quadrupole 
moment Qij of one star depends on the tidal field £y 
from the monopole of the other star through the rela- 
tion Qij — —X£ij, where A is the EOS dependent tidal 
deformability and is related to the neutron star's dimen- 
sionless Love number &2 and radius R through the re- 
lation A = g|jfc2i? 5 . The leading [l — 2) relativistic 
tidal Love number &2 was first calculated in Ref. [10] for 
polytropic EOS, then for EOS with hadronic and quark 



matter [TTJ [T2], as well as for EOS with analytic solu- 
tions to the stellar structure equations [12] ■ Its effect 
on the binary inspiral (including the contribution due 
to tidally excited f-modes) was calculated to leading or- 
der [13], and later extended to 1PN order [UCG3. The 
gravitoelectric and gravitomagnetic tidal Love numbers 
for higher multipoles were calculated in (THl \U\ ■ The en- 
ergy has now been calculated to 2PN order in the tidal 
corrections in the effective one body (EOB) formalism, 
including £ = 2 and 3 gravitoelectic interactions and the 
£ = 2 gravitomagnetic interaction, using the effective ac- 
tion approach [18], and most terms in the EOB wave- 
form are now known to 2.5PN order in the tidal interac- 
tions |19| . Finally, the accuracy of the adiabatic approxi- 
mation to tidal interactions was calculated using an affine 
model, and a Love function was found that corrects for 
this approximation and asymptotically approaches the 
Love number for large binary separations |201 121] . 

The measurability of tidal parameters by detectors 
with the sensitivity of aLIGO and ET was examined 
for BNS inspiral for gravitational wave frequencies below 
450Hz [T3] using polytropic EOS as well as for theoretical 
hadronic and quark matter EOS [IT] , The studies found 
that tidal interactions were observable during this early 
inspiral stage (prior to the last ~ 20 gravitational wave 
cycles before merger) only for stiff EOS and NS masses 
below 1.4 Mq. On the other hand, using tidal correc- 
tions up to 2.5PN order in the EOB approach, it was 
found that tidal parameters are in fact observable when 
including the extra ~ 20 gravitational wave cycles up to 
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the point of contact [T5] . 

Numerical simulations have also been used to study the 
measurability of matter effects during the late inspiral of 
BNS systems. Read et al. [55] examined the measurabil- 
ity of EOS information during the last few orbits, using 
numerical simulations assuming that non EOS parame- 
ters do not correlate significantly with EOS parameters. 
They found that the NS radius could be measured, us- 
ing only the last few orbits, to an accuracy of ~ 10%. 
In addition, using simulations in the conformally flat ap- 
proximation, Ref. [53] found that the NS radius could 
be determined to a comparable accuracy from the post- 
merger phase. 

Numerical work is also in progress to verify the ac- 
curacy of the inspiral tidal description using BNS sim- 
ulations. Initial studies comparing quasiequilibrium se- 
quences [53] and full hydrodynamic simulations [55] [55] 
with post-Newtonian and EOB tidal descriptions found 
noticeable differences, suggesting large corrections to the 
analytic description might be necessary. However, later 
comparisons have highlighted the importance of the hy- 
drodynamics treatment, numerical resolution, and wave- 
form extraction radius in determining the tidal contribu- 
tion to the waveform with numerical simulations [57] [55] . 
These later comparisons found that within numerical er- 
ror the improved simulations are consistent with the 2PN 
accurate EOB tidal description up to merger [27], or at 
least up to the last ~ 10 gravitational wave cycles [55] , 

Work is now underway to understand the measur- 
ability of EOS information in black hole-neutron star 
(BHNS) systems as well. Studies using the inspiral wave- 
form with tidal corrections up to 1PN order found that 
BHNS waveforms are not distinguishable from binary 
black hole (BBH) waveforms before the end of inspiral 
when the frequency reaches the ISCO frequency or tidal 
disruption frequency |29j . This conclusion was also found 
to be true for BHNS systems with spinning black holes 
where this ending frequency may be larger than for a 
nonspinning BH. On the other hand, work examining a 
possible cutoff in the gravitational wave amplitude due 
to tidal disruption of the neutron star by the black hole 
suggested that the NS radius may, in fact, be measurable 
with second generation detectors [30lI32| . 

Several numerical simulations have also been carried 
out to examine EOS dependent effects during the end of 
the BHNS inspiral [33- 36 , and these simulations now in- 
clude BH to NS mass ratios of up to Q = Mbh/M^s = 
7.1 and black hole spins of up to xbh = Jbh/Mbh = 0-9, 
where Jbh is the black hole's angular momentum |37U38| . 
In Ref. |39j (hereafter Paper I), we examined numerical 
simulations of the last few orbits, merger, and ringdown 
for systems with nonspinning black holes and low mass 
ratios of Q = 2 and 3. We found that when consider- 
ing only the merger and ringdown, the tidal deformabil- 
ity A was the best measured EOS parameter and was 
marginally measurable for second generation detectors. 
As in Paper I, we define a dimensionless version A of the 



the tidal deformability using the NS mass 



A := GA 



GM] 



NS 



~ 3 2 \GMns. 
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In this paper we repeat the analysis of Paper I for mass 
ratios up to Q — 5 and black hole spins from %bh = — 0.5 
to 0.75. We will also address many of the simplifications 
used in Paper I that can have a significant impact on the 
detectability of EOS parameters. Previously, we consid- 
ered only the tidal information that could be obtained 
from the merger and ringdown, ignoring the small accu- 
mulating phase drift during the inspiral that results from 
tidal interactions. We will find that coherently adding the 
slow tidal phase drift from the inspiral to the tidal effect 
during the merger and ringdown (as was also discussed 
in Ref. |38j ) can improve the measurability of tidal pa- 
rameters by a factor of ~ 3 over just the merger and 
ringdown. 

In Paper I we also ignored possible correlations be- 
tween the tidal parameter A and the other binary pa- 
rameters when estimating the measurability of A with the 
Fisher matrix approximation. We have addressed this is- 
sue by developing a frequency-domain, analytic BHNS 
waveform model fit to our BHNS simulations, and based 
on analytic BBH waveform models. This fit allows us to 
accurately evaluate derivatives in the Fisher matrix and 
evaluate correlations between the tidal parameter A and 
the other parameters. We find that although these cor- 
relations are not nearly as strong as with other pairs of 
parameters, they can increase the uncertainty in A by a 
factor of ~ 3. Overall, we find that the measurability of 
A using the improved methods in this paper are about 
the same as the estimates in Paper I, where only the 
merger and ringdown were considered and correlations 
with the other parameters were presumed to be negli- 
gible. Finally, we calculate the systematic error in our 
BHNS waveform model and find that it is significantly 
smaller than the statistical error for aLIGO and compa- 
rable to the statistical error for the ET detector. 

In Sec. [H] we briefly describe the EOS and numerical 
methods used to generate our BHNS waveforms, as well 
as the waveform behavior as a function of mass ratio, 
BH spin, and EOS. We then construct hybrid waveforms 
that join the numerical waveforms to inspiral waveforms 
in Sec. |III[ We then provide an overview of parameter 
estimation and show that A is the best measured EOS 
parameter in Sees. IV and[V] Finally, to understand cor- 
relations between the parameters, we construct a phe- 
nomenological BHNS waveform fit to our hybrid BHNS 



waveforms in Sec. VI then estimate the ability to mea- 



sure A with aLIGO and ET in Sec. IVHI We conclude in 
Sec. |VIlI] with a brief discussion of some of the improve- 
ments needed to generate templates accurate enough to 
use in data analysis pipelines. 

Conventions: We use the following sign convention for 
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the Fourier transform of a signal x(t) 

/•OO 

£(/)= / x(t)e +2 * ift dt, 



(2) 



and we will decompose the complex Fourier transform 
into amplitude and phase as h(f) = \h(f)\e + ^^\ These 
sign conventions are opposite those of Paper I, and are 
chosen to agree with those of the PhenomC waveform 
model [40 , which we use extensively to construct hybrid 
waveforms. In addition, we set G = c = 1 unless other- 
wise stated. 



II. SIMULATIONS 

Following Paper I on nonspinning BHNS systems, we 
perform a large set of simulations where we systemati- 
cally vary the parameters of a parametrized EOS, then 
look for the combination of parameters that are best 
extracted from gravitational-wave observations. Specifi- 
cally we choose a simplified two-parameter version of the 
piecewise polytrope EOS introduced in Ref. 41 . For 
this EOS, the pressure p in the rest-mass density interval 
Pi-i < p < pi is 



p(p) = K lP 



r, 



(3) 



where Ki is a constant, and Ti is the adiabatic index. We 
fix the crust EOS defined by densities below the transi- 
tion density p . In the crust, K = 3.5966 x 10 13 in 
cgs units and T = 1.3569, such that the pressure at 
10 13 g/cm 3 is 1.5689 x 10 31 dyne/cm 2 . Above the tran- 
sition density po, the core EOS is parametrized by the 
two parameters p\ and T\, The pressure p\ is defined as 
the pressure at pi — 10 14 7 g/cm 3 and the adiabatic in- 
dex Ti of the core will, for simplicity, be written T. The 
constant Ki for the core is then given by K\ = p\j ' p\. 
Once the two parameters p± and T are set, the divid- 
ing density pa between the crust and the core is given 
by the density where the crust and core EOS intersect: 
Pa = (AVAd) 1/(r ~ ro) . Finally, given this EOS, the en- 
ergy density e can be evaluated by integrating the first 
law of thermodynamics 



P 



-pd- 



,1 



(4) 



As discussed in more detail in Paper I, quasiequilib- 
rium configurations are used as initial data for the sim- 
ulations S2] j and are computed using the spectral- 
method library LORENE [33]. The numerical simula- 
tions are performed using the adaptive-mesh refinement 
code SACRA [44]. To obtain the gravitational wave- 
form h + — ih x , the outgoing part of the Weyl scalar 
^4 = h + — ih x is extracted from these simulations at a 
finite coordinate radius, and is then integrated twice us- 
ing a method known as Fixed Frequency Integration |45j . 
Specifically, we take the Fourier transform of ^4, then 



integrate twice in time by dividing by (27ri/) 2 . Low fre- 
quency components are filtered out as in Paper I, and the 
inverse Fourier transform is then taken to find h + — ih x 
in the time domain. 

We have performed 134 simulations of the late inspi- 
ral, merger, and ringdown of BHNS systems, using 21 
sets of parameters for our two-parameter EOS. We have 
also varied the mass ratio from Q = 2 to 5, the spin of 
the black hole from xbh = —0.5 to 0.75, and the neutron 
star mass from M NS = 1.20 M to 1.45 M . The EOS 
parameters used as well as the corresponding NS radius, 
Love number, and tidal deformability for the three NS 
masses can be found in Table [TJ A list of all the simula- 
tions and their starting frequencies is given in Table [il] In 
addition, we plot the EOS as points in parameter space in 
Fig. [T] along with contours of radius, tidal deformability 
A, and maximum NS mass. The 1.93 M maximum mass 
contour corresponds to the constraint from the recently 
observed pulsar with a mass of 1.97 ± 0.04 M measured 
using the Shapiro delay [35]. For this parametrized EOS, 
parameters below this curve have a maximum mass less 
than 1.93 M , and therefore do not agree with the mea- 
surement of a NS with M NS > 1.93 M Q . 



35.0 



34.8 



42.34.6^-" 



534.4 



34.2 



34.( 













R=12km 




A=500 


f \ 





2.5 



3.0 

r 



3.5 



4.0 



FIG. 1: The 21 EOS used in the simulations are repre- 
sented by blue points in the parameter space. For a NS of 
mass 1.35 M@, NS radius contours are solid blue and tidal 
deformability contours are dashed red. Also shown are two 
dotted contours of maximum NS mass. EOS parameters in 
the shaded region do not allow a 1.35 Mq NS. 

In Fig. [2] we show two representative waveforms. The 
waveform with a very soft EOS (p.3r2.4), and therefore 
small radius and tidal deformability, behaves very much 
like a BBH waveform where the inspiral smoothly tran- 
sitions to quasinormal mode ringdown. For the stiff EOS 
(p.7r3.0), however, the neutron star is tidally disrupted 
near the end of inspiral; the disruption and the spread of 
tidally stripped matter to form a roughly axisymmetric 
disk leads to a rapid decrease in the waveform amplitude 
and suppresses the subsequent ringdown. 
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TABLE I: Neutron star properties for the 21 EOS used in the simulations. The original EOS names 22 34, 36 are also listed, 
pi is given in units of dyne/cm 2 , maximum mass is in Mq, and neutron star radius R is in km. R, k2, and A are given for the 
three masses used: {1.20, 1.35, 1.45} Mq. 



EOS 


logpi 


r 


Mnax 


-Rl.20 


k 


2,1.20 


Al.20 


-Rl.35 


^2,1.35 


Al.35 


-Rl.45 


^2,1.45 


Ai.45 


p.3r2.4 


Bss 


34.3 


2 


4 


1.566 


10.66 





0765 


401 


10.27 


0.0585 


142 


9.89 


0.0455 


64 


p.3r2.7 


Bs 


34.3 


2 


7 


1.799 


10.88 





0910 


528 


10.74 


0.0751 


228 


10.61 


0.0645 


129 


P.3P3.0 


B 


34.3 


3 





2.002 


10.98 





1010 


614 


10.96 


0.0861 


288 


10.92 


0.0762 


176 


p.3r3.3 




34.3 


3 


3 


2.181 


11.04 





1083 


677 


11.09 


0.0941 


334 


11.10 


0.0847 


212 


pAT2A 


HBss 


34.4 


2 


4 


1.701 


11.74 





0886 


755 


11.45 


0.0723 


301 


11.19 


0.0610 


158 


pAT2.7 


HBs 


34.4 


2 


7 


1.925 


11.67 





1004 


828 


11.57 


0.0855 


375 


11.47 


0.0754 


222 


p.4r3.0 


HB 


34.4 


3 





2.122 


11.60 





1088 


872 


11.61 


0.0946 


422 


11.59 


0.0851 


263 


p.4r3.3 




34.4 


3 


3 


2.294 


11.55 





1151 


903 


11.62 


0.1013 


454 


11.65 


0.0921 


293 


p.5r2.4 




34.5 


2 


4 


1.848 


12.88 





1000 


1353 


12.64 


0.0850 


582 


12.44 


0.0747 


330 


p.5r2.7 




34.5 


2 


7 


2.061 


12.49 





1096 


1271 


12.42 


0.0954 


598 


12.35 


0.0859 


366 


p.5r3.0 


H 


34.5 


3 





2.249 


12.25 





1165 


1225 


12.27 


0.1029 


607 


12.27 


0.0937 


387 


p.5r3.3 




34.5 


3 


3 


2.413 


12.08 





1217 


1196 


12.17 


0.1085 


613 


12.21 


0.0995 


400 


p.6r2.4 




34.6 


2 


4 


2.007 


14.08 





1108 


2340 


13.89 


0.0970 


1061 


13.73 


0.0875 


633 


P.6T2.7 




34.6 


2 


7 


2.207 


13.35 





1184 


1920 


13.32 


0.1051 


932 


13.27 


0.0960 


585 


p.er3.o 




34.6 


3 





2.383 


12.92 





1240 


1704 


12.97 


0.1110 


862 


12.98 


0.1022 


558 


p.6r3.3 




34.6 


3 


3 


2.537 


12.63 





1282 


1575 


12.74 


0.1155 


819 


12.79 


0.1068 


541 


p.7r2.4 




34.7 


2 


4 


2.180 


15.35 





1210 


3941 


15.20 


0.1083 


1860 


15.07 


0.0995 


1147 


p.7r2.7 




34.7 


2 


7 


2.362 


14.26 





1269 


2859 


14.25 


0.1144 


1423 


14.22 


0.1058 


912 


p.7r3.0 


1.5H 


34.7 


3 





2.525 


13.62 





1313 


2351 


13.69 


0.1189 


1211 


13.72 


0.1104 


795 


p.7r3.3 




34.7 


3 


3 


2.669 


13.20 





1346 


2062 


13.32 


0.1223 


1087 


13.39 


0.1140 


726 


p.9r3.0 


2H 


34.9 


3 





2.834 


15.12 





1453 


4382 


15.22 


0.1342 


2324 


15.28 
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1560 



p.3r2.4 (R=10.27km, A=142) s\ 

: p.7r3.0 (R=13.69km, A=1211) ^^\\ 




l| fill Ilk 


A/WWW\j 

A-bh=0, Q=2, M NS =1 .35 M a , Mn =0.028 





200 400 600 800 

t/M 



FIG. 2: Polarization h+ evaluated on the orbital axis for 
two waveforms from numerical BHNS simulations that differ 
only in their EOS. Also shown in black are the waveform 
amplitudes \h+ ~ ih x \. 



Because trends in the BHNS waveform are most appar- 
ent in terms of the amplitude and phase of the Fourier 
transform, and because data analysis is usually done in 
the frequency domain, wc will now focus our discussion 
of the waveforms on the frequency-domain waveform be- 
havior. Several representative waveforms with varying 
tidal deformability A, mass ratio Q, and spin xbh are 
shown in Figs. [3}]7] 

As was found in Paper I, the waveform monotonically 
departs from a BBH (A = 0) waveform as A increases, 
and this is true for systems with spinning black holes as 



well, as we see from Figs. |3j|5| In particular, the cutoff 
frequency, where the waveform begins a sharp drop in 
the amplitude, decreases monotonically with increasing 
A. The accumulated BHNS phase <&bhns at fixed fre- 
quency / similarly decreases with increasing A, because 
the orbit loses energy more rapidly: There is less time for 
the phase to accumulate. As a result, the departure of 
^bhns from the accumulated BBH phase $bbh increases 
with increasing A. 

More massive black holes exert smaller tidal forces on 
their companion near coalescence, because the radius of 
the innermost orbit is roughly proportional to Mbh- As 
a result, the difference in amplitude and phase between a 
BHNS and BBH waveform decrease when the mass ratio 
Q increases. The effect is clear in Fig. [6j which displays 
the dramatically enhanced departure of amplitude and 
phase from that of a BBH waveform as Q decreases. On 
the other hand, the radius of the innermost orbit de- 
creases with increasing aligned BH spin xbh, implying a 
larger maximum tidal force for larger xbh- The resulting 
enhanced departure from a BBH waveform is shown in 
Fig.0 



5 



TABLE II: Data for the 134 BHNS simulations. NS mass is in units of Mq. Mfio is the initial orbital angular velocity of the 
system. The columns labeled "Phen" and "EOB" indicate if analytic PhenomC and EOB BBH waveforms are available for 
the given values of xbh and Q. "C" indicates the waveform is available and is calibrated for those parameter values. "NC" 
indicates the waveform could be generated but is not calibrated. "NA" indicates the waveform is not available. 



Xbh Q M ns 


EOS 


Mtt 


Phen EOB 


Xbh Q M ns 


EOS 


MQ 


Phen EOB 


Xbh Q M ns 


EOS 


MQ 


Phen EOB 


-0.5 


2 


1 


35 


p.3r3.0 


0.028 


C 


c 


0.25 


3 


1 


35 


p.7r3.0 


0.030 


C 


c 





75 


2 


1 


20 


P.3P3.0 


0.028 


C 


NA 


-0.5 


2 


1 


35 


p.4r3.0 


0.028 


C 


c 


0.25 


3 


1 


35 


p.9r3.0 


0.028 


C 


c 





75 


2 


1 


20 


p.4r3.0 


0.028 


C 


NA 


-0.5 


2 


1 


35 


p.5r3.0 


0.025 


C 


c 


0.25 


4 


1 


35 


p.3r3.0 


0.031 


C 


c 





75 


2 


1 


20 


p.5r3.0 


0.028 


C 


NA 


-0.5 


2 


1 


35 


p.9r3.0 


0.022 


C 


c 


0.25 


4 


1 


35 


p.5r3.0 


0.031 


C 


c 





75 


2 


1 


20 


p.9r3.0 


0.025 


C 


NA 


-0.5 


3 


1 


35 


p.4r3.0 


0.030 


C 


c 


0.25 


4 


1 


35 


p.7r3.0 


0.031 


C 


c 





75 


2 


1 


35 


p.3r3.0 


0.028 


C 


NA 





2 


1 


20 


p.3r3.0 


0.028 


C 


c 


0.25 


4 


1 


35 


p.9r3.0 


0.029 


C 


c 





75 


2 


1 


35 


p.4r3.0 


0.028 


C 


NA 





2 


1 


20 


p.4r3.0 


0.028 


C 


c 


0.25 


5 


1 


35 


p.3r3.0 


0.033 


NC 


c 





75 


2 


1 


35 


p.5r3.0 


0.028 


C 


NA 





2 


1 


20 


p.5r3.0 


0.028 


C 


c 


0.25 


5 


1 


35 


p.5r3.0 


0.033 


NC 


c 





75 


2 


1 


35 


p.7r3.0 


0.028 


C 


NA 





2 


1 


20 


p.9r3.0 


0.022 


C 


c 


0.25 


5 


1 


35 


p.7r3.0 


0.033 


NC 


c 





75 


2 


1 


35 


p.9r3.0 


0.025 


C 


NA 





2 


1 


35 


p.3r2.4 


0.028 


C 


c 


0.25 


5 


1 


35 


p.9r3.0 


0.033 


NC 


c 





75 


2 


1 


45 


p.3r3.0 


0.028 


C 


NA 





2 


1 


35 


p.3r2.7 


0.028 


C 


c 


0.5 


2 


1 


35 


p.3r3.0 


0.028 


C 


c 





75 


2 


1 


45 


p.4r3.0 


0.028 


C 


NA 





2 


1 


35 


p.3r3.0 


0.028 


c 


c 


0.5 


2 


1 


35 


p.4r3.0 


0.028 


c 


c 





75 


2 


1 


45 


p.5r3.0 


0.028 


C 


NA 





2 


1 


35 


p.3r3.3 


0.025 


c 


c 


0.5 


2 


1 


35 


p.5r3.0 


0.028 


c 


c 





75 


2 


1 


45 


p.9r3.0 


0.025 


C 


NA 
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FIG. 3: (Color online) Left panel: Amplitude \h\ of the Fourier transform of BHNS waveforms for xbh = 0.5, Q = 3, and 
Mns = 1.35Mq and 21 EOS. Color indicates the value of log(pi) while the line style indicates the value of V. Also shown are two 
analytic approximations to the BBH waveform with the same values of xbh and Q: PhenomC and EOB discussed in Section|TTl] 
Right panel: Phase of the Fourier transformed waveform relative to the PhenomC BBH waveform A<3> = $ — $phen. The curves 
are truncated when the amplitude D B g\h\/M drops below 0.1, and numerical error begins to dominate. The phase of the EOB 
BBH waveform is also shown relative to the PhenomC BBH waveform. In both figures the BHNS and EOB waveforms have 
been windowed, matched, and spliced to the tidally corrected PhenomC waveform as described in Section [Hi] The windowing 
width At w in = 300M; the start and end frequencies for matching, represented by solid vertical lines, are Mfi — 0.018 and 
Mff = 0.026; and the splicing interval, represented by dotted vertical lines, is Msi = Mfi (overlapping with the solid line) 
and Ms f = Ms; + 0.001. 
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FIG. 4: Same as Fig. [3] but X bh = -0.5, Q = 2, and M NS = 1.35M e 
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FIG. 5: Same as Fig. [3] but xbh = 0.5, Q = 4, and Mns = 1.35Mq. During the ringdown, the PhenomC waveform amplitude 
is noticeably less than the BHNS waveform with the softest EOS (p.3P3.0) around Mf = 0.1. The EOB amplitude, however, 
is always greater than the BHNS amplitude. 
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FIG. 6: BHNS waveforms with mass ratios of Q = {2,3,4,5} with all other parameters fixed at ^bh = 0, Mns == 1.35M Q , 
and the stiffest EOS = p.9r3.0. Left panel: Amplitude of BHNS waveform as well as PhenomC BBH waveforms (black dots) 
with the same values of Q and Xbh- The difference in amplitude between BBH and BHNS waveform decreases with the mass 
ratio Q. Right panel: Difference in phase between BHNS and PhenomC BBH waveform also decreases with Q. Windowing, 
matching, and splicing are identical to Fig. [3] 
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FIG. 7: Same as Fig. [6]except BH spin is varied instead of mass ratio. Waveforms with spins of xbh = {—0.5, 0, 0.5, 0.75} are 
shown, and Q = 2, A/ns = 1.35MQ, EOS = p.9r3.0 for all waveforms. In contrast with the mass ratio (Fig. 16b, the departure 
from the BBH waveform increases as the spin increases. However, the effect is moderate. 



III. CONSTRUCTING HYBRID 
INSPIRAL-MERGER-RINGDOWN WAVEFORMS 

To obtain as much information as possible about the 
physical parameters of a BHNS coalescence, we will con- 
struct a hybrid waveform that joins the numerical merger 
and ringdown waveforms to an analytic BHNS inspiral 
waveform. This inspiral waveform needs to account for 
the aligned spin xbh of the BH as well as tidal interac- 
tions through the parameter A. In addition, because the 
numerical waveforms only include the last ~ 10 GW cy- 
cles before merger, the inspiral waveform, including spin 
and tidal corrections, should be valid as close to merger as 
possible. Finally, when we construct an analytic BHNS 
inspiral-merger-ringdown (IMR) waveform model in Sec- 
tion |VI| we will find it useful to proceed by modifying 
a BBH IMR waveform, and so the waveform model will 
need to accurately model the merger and ringdown of a 
BBH system as well. 

We will use two classes of waveform models that have 
been calibrated to numerical BBH simulations. The pri- 
mary waveform model, described in the next subsection, 
is the frequency- domain BBH waveform labeled Phe- 
nomC [40], and we will add a 1PN tidal correction to 
the inspiral portion of this waveform. We will also use 
a time-domain effective one body (EOB) waveform that 
incorporates spin corrections [27], and we will again add 
tidal corrections to the inspiral. We compare these two 
waveform models with no inspiral tidal correction to each 
other in Figs.|3}|5j We note that although these two wave- 
forms agree well for small mass ratios and spin, their dif- 
ferences become important as the mass ratio and spin 
increase. 



A. PhenomC waveform with tidal corrections 

Several frequency-domain phenomcnological models 
are now available for the complete IMR BBH waveform. 
These models include the PhenomA [35] model for non- 
spinning BBH systems, as well as the PhenomB [49] 
model and improved PhenomC [30] model for aligned- 
spin BBH systems which we will use. 

In the PhenomC waveform, the Fourier transform 
of the waveform is decomposed into an amplitude 
Aphcn(Mf) and phase $ ph en(^/) as 



^phen 



(Mf) = A phcn (Mf)e 



AMf) 



(5) 



The inspiral is described by the stationary phase approxi- 
mation TaylorF2 post-Newtonian waveform, and the spin 
contribution to the waveform of both bodies is approxi- 
mated by the mass-weighted average spin [40] 



Xavg 



Mi M 2 



(6) 



For our BHNS systems where the NS spin is assumed 
to be negligible compared to the BH spin, x avg = 



(Msb./M)xbh- The amplitude and phase of the merger 
and ringdown for frequencies above Mf = 0.01 are then 
fit to a large set of numerical simulations with mass 
ratios from 1-4 and various combinations of aligned or 
anti-aligned BH spins up to \xi\ = 0.85 as described in 
Ref. [40]. 

We add tidal corrections to the inspiral for this wave- 
form using the TaylorF2 stationary phase approximation 
up to 1PN order [15]. This is the same quantity used by 
Pannarale et al. [H] who found that BBH and BHNS 
waveforms are indistinguishable by aLIGO when only 
considering the inspiral. Explicitly, we add a tidal cor- 
rection term iPt{M f) to the BBH phase, 

MMf) = t4^ M/) ~ 5/3 [ao(^M/) 10/3 + ai(^M/) 12/3 

(7) 

where in terms of the symmetric mass ratio 77 = 
M B hM ns /(Mbh + Mns) 2 and for M BH > M NS , 



a 
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(8) 



A. 



hi. 



insp(M/) = 
phen+T(Mf) = 



inspiral waveform is now 
A>hon(M/)e l *p h -+ TW \ where $ 

% h en(Mf)+MMf). 

We will work mostly with this frequency-domain wave- 
form in the sections that follow. Since a BHNS waveform 
enters the detector band starting at frequencies as low as 
10 Hz for aLIGO and 1 Hz for ET, it is much more effi- 
cient to start with a frequency-domain inspiral waveform 
than to evaluate the Fourier transform of a time-domain 
waveform. 



B. Spinning EOB waveform with tidal corrections 

The other waveform model we will use, time-domain 
EOB waveforms, have proven succesful at reproducing 
the complete IMR of nonspinning BBH waveforms [50] . 
and we used these EOB waveforms in Paper I to generate 
hybrid waveforms that did not incorporate tidal correc- 
tions for the inspiral. (See Appendix C of Paper I and 
references therein for a review of the EOB formalism.) 
Recently, spin terms have been calculated for the EOB 
Hamiltonian and resummed waveforms, and free param- 
eters for the merger have been calibrated to numerical 
nonspinning BBH waveforms for mass ratios from 1-6, as 
well as for equal mass, aligned-spin waveforms with spins 
of xi — Xi — ±0.44 [47] . We will use EOB waveform ta- 
bles generated by Taracchini and Buonanno [51 j - These 
tables were generated by evolving the EOB equations of 
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motion with an initial radial coordinate of r = 40M and 
a value of the radial velocity r consistent with the radi- 
ation reaction force to minimize initial eccentricity [52"] . 
The waveform is then evaluated starting at r = 30M 
where any residual eccentricity is negligible |51) . 

In addition, tidal interactions have also recently been 
incorporated in the EOB formalism. In the method pro- 
posed in Ref. [19] . a term representing the conservative 
part of the tidal interaction is added to the radial po- 
tential A(r) in the EOB Hamiltonian. Tidal corrections 
are also added to the resummed waveform hg m which are 
used to calculate the radiation reaction force in the equa- 
tions of motion. The solutions to the equations of motion 
are then plugged back into hg m to produce a final wave- 
form as a function of time. However, for simplicity and 
because the versions of the EOB formalism that incorpo- 
rate spin and tidal interactions are slightly different and 
have not been calibrated to simulations with both spin 
and matter, we will instead incorporate tidal interactions 
in the spinning EOB waveform using the same method as 
that for the PhenomC waveform. Specifically, we Fourier 
transform the EOB waveform, decompose it into ampli- 
tude -Aeob(-W/) and phase $eob(-W/), and then simply 
add the expression ifa(Mf) from Eq. Q to $ EOB (M/). 
Future work on spinning BHNS systems should treat the 
spin and tidal interactions consistently. 



and we choose the start of the windowing to be the start 
of the numerical waveform at w; = 0. This window- 
ing minimizes the oscillatory Gibbs phenomenon that re- 
sults from Fourier transforming a waveform segment with 
nonzero starting amplitude. When matching waveforms, 
a time constant r and phase constant 4> are free param- 
eters that need to be fixed. For a generic waveform h(t), 
the time and phase can be adjusted to produce a shifted 
waveform h shlit (t; r, <f>) — h(t — r)e 1 ^. The Fourier trans- 
formed waveform, which can be written in terms of ampli- 
tude and phase as h(f) = \h{f)\e l ^^\ has a correspond- 
ing shifted waveform /i shift (/;r,0) = \h{f)\e i9 ^^' T ^ , 
where 4> shift (/;r,0) = $(/)+2tt/t+0. When joining the 
inspiral and numerical waveforms, we hold the phase of 
the tidally corrected inspiral waveform < I > bbh-i-t(/) fixed 
and adjust the phase of the numerical waveform $nr(/), 
such that $fjf*(/; r, 0) = $nr(/) + 2vr/r + <f>. We then 
match the waveforms by performing a least-squares fit in 
the matching interval fi < f < ff (width A/ matc h) that 
minimizes the quantity 



J Vm (/;T,^)-$BBH+T(/)] 2 rff (10) 



C. Hybridization in frequency domain 

In Paper I, where we examined EOS effects only during 
the BHNS merger and ringdown, we used a time-domain 
matching method. In this paper, however, because we 
examine tidal effects for the entire IMR waveform, we 
will find it convenient to start with the inspiral waveform 
in the frequency domain before matching. The method 
we use closely follows the frequency-domain least-squares 
method used in Ref. [?D] for BBH waveforms. 

We begin by windowing the numerical BHNS waveform 
with a Hann window over the interval Wj to w/ (width 



Wf - Wi) 



w on (t) = 



1 



1 — cos 



ir[t- 

Wf 



Wi. 
- Wi 



to determine the free parameters r and (f>. 

Once the time and phase shifts are found, we smoothly 
turn on the numerical waveform and smoothly turn off 
the phenomenological waveform within a splicing window 
Si < f < s f (width A/ S pii ce ) using Hann windows 



Wos(f) 

Won(f) 



1 + COS 

1 — cos 



7T[/ 



Sf - Si 

7r[/ - Si] 

Sf - Si 



(11) 

(12) 



(9) The amplitude of the hybrid waveform is then 



and the phase is 



|^bbh(/)|, f<Si 

'l.vl,rid(/)| = { W oS (f)\h B Bu(f)\ +«W(/)|/INr(/)|, Si<f<Sf (13) 

I^nr(/)|, f>Sf, 



$BBH+t(/), / < Si 

$hybrid(/)= \ W off (/)$BBH+T(/) + Won(/)[ < 5NR(/)+2^/T + </>], Sl < f < S f (14) 

$nr(/) + 2tt/t + 4>, f>s f . 
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A hybrid waveform for the system (xbh = 0,Q = 
2,M NS = 1.35 M 0) EOS = p.5r3.0) is shown in Fig. § 
where we matched the numerical waveform to the Phe- 
nomC waveform with and without the inspiral tidal cor- 
rection ipr- In the right panel of Fig. [5] we show four 
phases relative to the PhenomC BBH phase $bbh. The 
black dashed curve is the PhenomC BBH phase. The 
solid black curve is the numerical BHNS waveform phase 
after it is matched directly to the BBH waveform as was 
done in Paper I 1 . The red dashed curve represents the 
phase of the tidally corrected inspiral waveform 'I'bbh+t, 
and the curve is given by the analytic expression for the 
inspiral tidal correction $bbh+t(A^/) — ®bbh(M f) = 
ipr(Mf). Finally, the solid red curve is the numeri- 
cal BHNS waveform matched to the tidally corrected 
inspiral waveform. We note that the difference after 
the matching window between the solid red and solid 
black curves depends only on the difference in the match- 
ing term 2tt(M /)(t /M) + <f> in the hybrid phase ^hybrid 
(Eq. (14)) between when the inspiral tidal correction 
ipT is (red curve) or is not (black curve) included. Be- 
cause this term is linear, the difference between the two 
curves grows linearly. This difference is approximately 
V , T(A^/mid) + (M/-Af/ mid )^(M/ mi d), where M/ mid is 
the midpoint of the matching interval and the ' indicates 
a derivative with respect to Mf. As we will see below, 
this linear term has a large impact on the measurability 
of tidal parameters. 



D. Sensitivity of hybrid waveform to matching 
parameters 

In the hybridization procedure described above, we are 
free to choose the window width Ai w j n used in Eq. Q as 
well as the matching interval AM/ matc h = Mff — M fi 
and midpoint of the matching interval M/ m jd = (M fi + 
Mff)/2 used in Eq. ( 10 ). If the numerical waveform were 
long and identical to the inspiral waveform within some 
interval, the choice of these free parameters would have 
no impact on the values of r and (/). However, there are 
several sources of error. Because the numerical waveform 
has finite length, the beginning of the waveform needs to 
be windowed before Fourier transforming the waveform 
to reduce the Gibbs phenomenon. The matching interval 
should exclude as much of the beginning of the waveform 
as possible because the numerical simulation takes time 
to settle down from inexact initial conditions which in- 
cludes some initial eccentricity. It should also exclude the 
merger and ringdown which cannot be described by the 
tidal terms for the analytic inspiral waveform. On the 
other hand, the matching window must be wide enough 



to average over ringing from the Gibbs phenomenon that 
remains after windowing, the effects of eccentricity in the 
simulations, and other numerical noise. To isolate the ef- 
fect that each of these waveform errors has on the values 
of t and <p, we will introduce them sequentially to the 
waveform that is matched to the inspiral waveform. 

We first examine the effect of the Gibbs phenomenon, 
present in finite length waveforms, on the time and phase 
shifts r and <fi. To do this, we begin with an (effectively) 
infinite length time-domain EOB waveform that includes 
the inspiral, merger, and ringdown. We then mimic a 
numerical waveform by making a truncated copy of this 
EOB waveform that starts ~ 10 GW cycles (800M) be- 
fore merger. We window the first Ai w ; n and Fourier 
transform this truncated waveform, then match it to the 
Fourier transformed original waveform such that the only 
matching error is due to the Gibbs phenomenon. In the 
top panel of Fig. [9] we see the post- matching hybrid phase 
at M f = 0.05 depends on the window width Ai win and 
the frequency interval defined by AM/ matc h and M/ m jd. 
However, if we increase either the window width (from 
At w i n = lOOil/ to 300Af) or match over a larger fre- 
quency interval (from AA// matcll = 0.002 to 0.008), we 
can reduce the dependence of the hybrid phase on the 
Gibbs Phenomenon. 

We next consider the effect of eccentricity on the hy- 
brid phase in the bottom panel of Fig. [9] The numerical 
BHNS simulations begin with quasicircular (zero radial 
velocity) initial conditions that ignore the small radial 
velocity due to radiation reaction. As a result, the inex- 
act initial conditions lead to a small initial eccentricity 
(eo ~ 0.03), which eventually dies down after several or- 
bits. We can mimick this effect by generating EOB wave- 
forms with equivalent eccentricity by starting the EOB 
equations of motion with the same quasicircular (zero ra- 
dial velocity) initial conditions as the simulations found 
by ignoring the radiation reaction term in the EOB equa- 
tions of motion 52 . We match an EOB waveform with 
the quasicircular initial conditions Mfio = 0.028 to an 
effectively infinite length, zero eccentricity EOB wave- 
form with otherwise identical parameters. As with the 
top panel, the eccentric EOB waveform exhibits Gibbs 
oscillations because it has a finite length, and this effect 
can be reduced by increasing the window width and fre- 
quency matching interval. There is also an additional 
offset that results from the initial eccentricity, and this 
offset eventually dies down around A// m ;d ~ 0.03. 



1 Paper I used a time-domain method to match a numerical BHNS 
waveform to an EOB BBH inspiral waveform, and the time- 
domain matching interval used in Paper I corresponds to a fre- 
quency interval slightly less than that shown here. 



Finally, in Fig. 10 we match two numerical BHNS sim- 
ulations to PhenomC inspiral waveforms with tidal cor- 
rections to generate a full BHNS IMR hybrid. If the 
inspiral waveform and numerical waveform are identical 
within a frequency interval, then the hybrid phase will 
be independent of the matching region within that inter- 
val, and there will therefore be a plateau in the curve in 
Fig. 10 We therefore identify the best matching region 
as the region centered on the maxima M/ m ;d ~ 0.022 in 
the top panel and M/ m ;d ~ 0.020 in the bottom panel. 
In addition, as in Fig. [9j increasing the window width 
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FIG. 8: Amplitude D eB \h(Mf)\/M (left) and phase <E>(M/) (right) for a numerical BHNS waveform matched to the PhenomC 
BBH waveform with and without the tidal correction ipT (Eq. Q). The waveform parameters are (xbh = 0, Q = 2, Mns = 
1.35 M0,EOS = p.5r3.0). The matching window /» < / < // is bounded by solid vertical lines, and the splicing window 
Si < / < Sf, which begins at s, = /j, is bounded by dotted vertical lines. Note that matching the numerical BHNS waveform to 
a BBH waveform without tidal corrections, as was done in Paper I, results in ignoring a phase term that accumulates linearly 
even after the matching region, and underestimates the effect of matter. We truncate the waveform when the amplitude drops 
below 0.1 denoted by a black dot. 



and matching interval reduces oscillations due to Gibbs 
phenomena. 

over- 



We note that the best matching region in Fig. 10 



laps somewhat with frequencies where eccentricity still 
effects the hybrid phase as seen in Fig. [9j We can move 
the matching region to slightly higher frequencies; how- 
ever, the analytic tidal correction will rapidly become 
inaccurate. We can estimate roughly the error of the in- 
spiral tidal phase term in the bottom panel of Fig. [10] as 
the ratio of the 1PN to leading tidal corrections (defined 
by iPipn/iPopn = (ai/a )(nMf) 2 / 3 ). 

Although the best choice of matching frequencies vary 
slightly depending on xbh, Q, and the EOS, setting the 
matching region to be the same for all waveforms changes 
the final hybrid phase no more than the matching uncer- 
tainties described above. From the above considerations, 
we choose as our windowing width At W i n = 300M, and 
matching parameters AM/ matc h = 0.008 and M/ mi d = 
0.022 such that the start and end matching frequencies 
are Mf z = 0.018 and Mf f = 0.026. In the matching re- 
gion, the amplitude and phase of the inspiral and shifted 
numerical waveforms agree reasonably well, so the choice 
of Si and sj does not significantly effect the results. We 
choose Msi = Mf l and Ms/ = Msi + 0.001. 



IV. PARAMETER ESTIMATION 

The primary goal of this paper is to determine how ac- 
curately EOS parameters can be measured from BHNS 
observations. In this section we will discuss the statis- 
tical error associated with detector noise as well as the 
systematic error that results from using inexact waveform 
templates to estimate parameters. In the next three sec- 
tions, we will then determine the combination of EOS pa- 
rameters that is best measured and generate an analytic 



waveform model matched to numerical waveforms that 
we then use to estimate the statistical and systematic 
errors in measuring that EOS parameter combination. 

The output of a gravitational-wave detector s(t; 6?) 
ri: i . + ft, B (t;0 T ) is the sum of detector noise n(t) and 
a possible gravitational-wave signal exactly described by 
fiE(t;0T) with the true parameters 9t- We assume the 
noise is a stationary, Gaussian time series, and therefore 
characterized by its power spectral density (PSD) 5 n (|/|) 
defined by the ensemble average 



1 



(n(f)n*(f)) = ^(f-f)S n (\f\), 



and its probability distribution 



p n [n(t)] oc 



,-{n,n)/2 



(15) 



(16) 



Here, (a, b) is the usual inner product between two time 
series a(t) and b(t) weighted by the PSD 

M = iSe J-Wl 4 . (17) 

The gravitational wave signal is given in terms of the two 
polarizations of the gravitational wave by 



h E (t; 0t) = F+h E+ (t; T ) + F x h Ex [t; T ), 



(18) 



where i*+,x are the detector response functions and de- 
pend on the location of the binary and the polarization 
angle of the waves. As in Paper I, we assume the binary 
is optimally located at the zenith of the detector and 
optimally oriented with its orbital axis along the line of 
sight. This condition is equivalent to averaging h + and 
h x (F + = F x = 1/2). 

In searches for gravitational- wave signals from compact 
binary mergers, a set of templates hit; 9) with parameters 
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FIG. 9: Dependence of hybrid waveform phase at Mf = 0.05 
on window width Af wm /M and matching interval with width 
A M /match and midpoint M/ m id. Top panel: An EOB wave- 
form, truncated to only include the last 800M before merger, 
is Fourier transformed then matched to the same EOB wave- 
form that has not been truncated. The oscillations in the hy- 
brid phase result from the Gibbs phenomenon. This can be 
partially suppressed by increasing the window width Ai wm /M 
as well as the width of the matching window AM/„ k h. Bot- 
tom panel: The effect of eccentricity can be evaluated by 
matching an eccentric EOB waveform to a long, zero ec- 
centricity EOB waveform. Here, an eccentric EOB wave- 
form is generated by using the quasicircular initial conditions 
MQo = 0.028. The dependence of the hybrid phase on the 
matching interval can be reduced by using a larger match- 
ing window AM/ mllc h. For both panels, xbh = and Q — 2. 
Because the overall phase of the hybrid waveform is arbitrary, 
we have set it to in this figure when Mf = 0.05. 



9 are compared to the signal s(t). The parameters that 
maximize the signal to noise ratio (SNR) 



(M) 



(19) 
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FIG. 10: Dependence of hybrid waveform phase at Mf = 
0.05 on window width Ai w i n /M and matching interval as 
in Fig. [9] Top panel: BHNS simulation with parameters 
{Xbh = 0, Q = 2, M NS = 1.35Af , EOS = p.3r2.4} is 
matched to the PhenomC BBH waveform with tidal phase 
corrections. Middle panel: BHNS simulation with parameters 
{xbh = 0, Q = 2, Mm = 1.35M Q , EOS = p.9r3.0}. Bottom 
panel: Relative contribution of the 1PN tidal phase correction 
to the leading order tidal phase correction provides a crude 
estimate of the error in the tidal phase correction. For refer- 
ence, dashed vertical lines give the ISCO frequency l/(6 3 ^ 2 7r) 
as well as the innermost circular orbit for nonrotating black 
holes, defined by the minimum of the 3PN energy (Eq. (194) 
of [53]), for mass ratios of Q = 2 and 5. 



are the best estimate of the true parameters 9t- For a 
template hE(t',9) that exactly represents the true wave- 
form, we will denote the best estimate of the true param- 
eters produced by the exact template by 6e- In practice, 
however, we only have an approximate template ft^i; 9), 
and we will denote the best estimate produced by this 
approximate template by 9a- 



In the large SNR limit, the difference Ad = 6a — 9t be- 
tween the best estimate using an approximate template 
and the true parameters of the binary system obeys a 
Gaussian distribution [53]. Specifically, for N parame- 
ters, the conditional probability of the error A9 given 
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the best estimate 9 A is 



V(2^) 7V det(E„) 



e -is-. 1 (A0 I -<Ae I ))(Ae J -<Ae J ))^ 



The mean is found to be approximately 



where 



(AO,) ^ -T~ i 1 (Sh(e A ),d ] h A (9 A )), 



Tij = {d l h A (9 A ),d J h A {9A)) 



(20) 



(21) 



(22) 



is the Fisher matrix and Sh(t; 9) — h A (t; 9) — 9) is 
the difference between the approximate and exact wave- 
form templates. The covariance between the parameters 
is [54] 

Etf = ((A0i - (A^»(A^ - (A^-») = rr 1 , (23) 

and the variance in A9i is therefore 

aN((A^-(A^)) 2 )=r^ 1 , (24) 

where the repeated indices in T^ 1 do not represent sum- 
mation. The statistical error ellipsoid to n standard de- 
viations is a contour of p(A9\9 A ) given by 



(A9 i -(A9 i ))(A9 j -(A9 j ))^ = 



(25) 



In addition, we identify the quantity A9 ayst = (A9) as 
the systematic error that results from using an approxi- 
mate waveform template instead of the exact waveform 
template. We will use this expression below as a criteria 
for the accuracy of our analytic BHNS waveform model. 



V. BEST MEASURED EOS PARAMETER 

In Paper I we found that, during the merger and ring- 
down, the best-measured combination of EOS parame- 
ters for nonspinning BHNS systems was consistent with 
the tidal deformability A. We used hybrid waveforms 
that ignored the inspiral tidal correction i/jt (Eq. d7|)), 
and only included EOS information from the merger and 
ringdown of the numerical part of the waveform. Wc 
then evaluated a restricted two-parameter Fisher matrix 
for the EOS parameters log(pi) and T, ignoring possible 
correlations between the EOS parameters and the other 
parameters. In this section we compare the measurabil- 
ity of EOS parameters for only the merger and ringdown 
to a waveform that includes EOS information in the full 
IMR hybrid waveform, and we do this for three combi- 
nations of mass ratio and black hole spin. In the next 
two sections we will address the issue of correlations be- 
tween EOS and non-EOS parameters by constructing an 
analytic BHNS waveform and calculating the complete 
Fisher matrix. 

As in Paper I, we evaluate the Fisher matrix from a 
set of hybrid waveforms by differentiating the waveform 



with respect to each parameter using finite differencing 
with two or more waveforms for each parameter. We fol- 
low the third method in Appendix A of Paper I which 
results in the greatest accuracy given the phase differ- 
ence between waveforms which can be several radians 
for EOS parameters. Specifically, we decompose each 
Fourier transformed hybrid waveform into the log of the 
amplitude In A(f;9) and accumulated phase &(f;9) 

/.(/j^e 11 *^, (26) 

then evaluate di In A and 9j$ individually. The derivative 
is now approximated by 



dih(f;6)tte 



inA(m+iHm I AlnA (f;0) , - A$(/; 



A9 l 



A9 l 



(27) 

where A/A9 1 represents central differencing, and In A 
and $ are evaluated at the midpoint with linear interpo- 
lation. 

Calculating the complete Fisher matrix using hybrid 
waveforms requires one to evaluate partial derivatives 
with respect to all parameters at a single point. For an 
aligned-spin BHNS system with 2 EOS parameters and 
a single detector, the waveform will have the form 

Msff (28) 

x e i[2wft c +^ c +g^(f;M,r,, X BH,log( Pl ),r)] 

where g A and g$, are generic functions, and there are 
8 parameters. The 5 intrinsic parameters are the chirp 
mass M. — (A/bh-Mns) 3 ^ 5 /-^ 1 ^ 5 : symmetric mass ratio 
T) = MbhM N s/M 2 , black hole spin xbh 2 , and the 2 EOS 
parameters log(pi) and T. The 3 extrinsic parameters, 
which can be differentiated analytically, are time of coa- 
lescence t c , phase of coalescence (j> c , and an effective dis- 
tance D c g that incorporates the true distance D as well 
as the orientation and sky location of the binary. (For 
an optimally oriented and located binary, D e f{ = D). If 
using central differencing, this requires 10 hybrid wave- 
forms for the 5 numerical derivatives at each point in 
the waveform parameter space, and is computationally 
expensive if one wants to explore the entire parameter 
space. In addition, in contrast to the small EOS depen- 
dent effects, small changes in Ai, T), and xbh can result in 



2 There will also be a neutron-star spin contribution XNS- Mag- 
netic dipole radiation, however, is expected to spin down a NS 
to a small fraction of the Kepler frequency well before the binary 
reaches the detector band. Furthermore, for the PhenomC wave- 
form, the aligned spins of the two bodies are approximated by 
the single mass- weighted average-spin parameter Xavg (Eq. |6]l). 
Since we will use Xbh as our onr y spin parameter it effectively 
becomes a linear combination of the BH and NS spins. Incorpo- 
rating NS with significant spins would require a different inspiral 
waveform model with two separate spin parameters as well as 
BHNS simulations with spining NS. 
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a large change in the phase of the waveform. This means 
that the simulations must be closely spaced in parameter 
space in order to accurately calculate derivatives, requir- 
ing a very large number of waveforms. In this section, 
we will therefore restrict the Fisher matrix calculation to 
the two EOS parameters log(pi) and T, and will leave to 
the next section a better way to differentiate the other 
parameters. 

For the BHNS systems discussed here, the greatest de- 
parture from BBH behavior occurs for gravitational- wave 
frequencies in the range 300-3000 Hz. As a result, de- 
tector configurations optimized for detection of BHNS 
systems with low noise in the region below 300 Hz may 
not optimally estimate EOS parameters. We therefore 
present results for the broadband aLIGO noise curve [57J 
and the ET-D noise curve [55] shown in Fig. 11 The 



shown in Fig. 

broadband aLIGO configuration uses zero-detuning of 
the signal recycling mirror and a high laser power, re- 
sulting in significantly lower noise above 300 Hz at the 
expense of slightly higher noise at lower frequencies. Sev- 
eral configurations have been considered for the Einstein 
Telescope denoted ET-B [55], ET-C [BUJ, and ET-D [55]. 
We use the most recent ET-D configuration and note that 
in the 300-3000 Hz range all of the ET configurations 
have a similar sensitivity. The published noise curves, 
and those used in this paper, are for a single interferom- 
eter of 10 km with a 90° opening angle. The current ET 
proposal is to have three individual interferometers each 
with a 60° opening angle configured in an equilateral tri- 
angle. This will shift the noise curve down appoximately 
20% [551. 



Q=2, xbh=0, A=0 

— Q=2, xbh=0, A=2000 

Q=5, ^bh=0.75, A=0 

— - Q=5, ;tbh=0.75, A=2000 
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FIG. 11: (Color online) Noise PSD for broadband aLIGO, 
ET-B, and ET-D. Also shown are the weighted amplitudes 
2f 1 / 2 \h\ of phenomenological waveforms for 4 parameter val- 
ues. For all waveforms, Mns = 1.35M0 and D e g = 100 Mpc. 
Black circles represent the start of the waveform fit at Mf = 
0.01. 

The l-<7 uncertainty ellipses in the EOS parameter 
space 6 = {log(pi),L} are A^A^Ty = 1. When cal- 
culating these 2-parameter error ellipses with the Fisher 



matrix using finite differencing, we sometimes find that 
two waveforms with the same log(pi) but different L are 
nearly identical, leading to derivatives that are domi- 
nated by numerical errors in the waveforms. To avoid 
this problem, we transform the L-log(pi) coordinate sys- 
tem to the u-v coordinate system shown in Fig. |12| We 
then evaluate the derivatives in this new coordinate sys- 
tem, and finally transform back with the chain rule 
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FIG. 12: Example of coordinates in the u-v coordinate sys- 
tem used to compute derivatives in Eqs. ( |29[ ) and ( |30[ |. The 
u-v coordinate system is chosen so the axes are not aligned 
with contours of A 1 / 5 . 



These ellipses are shown in Fig. 13 for the ET-D noise 
PSD for both nonspinning and spinning simulations when 
the BHNS waveform is matched to a PhenomC BBH 
inspiral waveform with no tidal correction ifo. As in 
Paper I we find the uncertainty contours are approxi- 
mately aligned with tidal deformability contours A 1//5 , 
and this holds for systems with spinning black holes as 
well. As in Paper I we plot A 1 / 5 instead of A because it 
is more closely related to the NS radius. We also note 
that the error ellipses found here for (xbh = 0, Q = 2, 
-^ns = 1.35Mq) using the PhenomC BBH inspiral wave- 
form and frequency-domain match are very similar to the 
results found in Fig. 11 of Paper I where we used the EOB 
BBH inspiral waveform and a time-domain match. 

In contrast, when the tidal correction ipj- is added to 
the PhenomC inspiral waveform before generating a hy- 
brid waveform, there is an improvement of roughly a fac- 
tor of 3 in the measurability of A as shown in Fig.[l4] The 
majority of the improvement arises because, as stated in 
the discussion of Fig. |8j even though the inspiral tidal 
correction is small, the hybridization procedure also adds 
a tidal term that grows linearly with frequency to the 
merger and ringdown which is not present when the nu- 
merical waveform is joined to an inspiral waveform with- 
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FIG. 13: l-cr error ellipses for the 2-parameter Fisher ma- 
trix using the ET-D noise curve. We use hybrid waveforms 
where the numerical BHNS waveform is matched to a Phe- 
nomC BBH inspiral waveform with no tidal correction. Bi- 
nary is optimally oriented and at a distance of 100 Mpc. Val- 
ues of Xbh, Q, and Mns are listed in each panel. Evenly 
spaced contours of constant A 1 ' 6 are also shown. Each ellipse 
is centered on the estimated parameter 6a denoted by a x . 
Top: Matching window has width AM/ matc h = 0.008 and is 
centered on M/ mid = 0.016. Middle: AM/ matc h = 0.008 
and M/ mid = 0.016. Bottom: AM/ match = 0.008 and 
M/ mid = 0.020. 



out tidal corrections. In addition, because the inspiral 
tidal correction ipx an d the tidal contribution to the 



matching term [fa (M f mid ) + (Mf-M / mid ) V4 ( M /mid )] 
are analytically proportional to A, the ellipses align much 
more closely with the A contours. We emphasize that the 
majority of the improvement comes from the above tidal 
contribution to the matching term, and not from the in- 
spiral term fa. As already noted [29], tidal interactions 
during the inspiral alone are not separately measurable. 
Finally we note that when using the broadband aLIGO 
noise curve instead of the ET-D noise curve, the error el- 
lipses in Figs. |l~3"| and [14| have nearly identical shape and 
orientation but the size is a factor of ~ 10 larger. 

We tentatively conclude that A is the dominant EOS 
dependent quantity that can be measured for the merger 
and ringdown as well as for the inspiral. We have consid- 
ered only the two-dimensional cross section Adjjjgg = of 
the eight-dimensional ellipsoid A^A^T.y = 1, however, 
and have therefore ignored correlations between the EOS 
parameters and the non-EOS parameters. We partially 
address this deficiency in the next section by explicitly 
accounting for correlations between A and the non-EOS 
parameters on a slice of constant T. One may still worry 
about the effect of correlations between T and all other 
parameters. Because the majority of the tidal contri- 
bution comes from the linearly growing matching term 
that analytically depends on A, we believe this assump- 
tion is mostly justified. We will test this assumption by 
calculating the systematic error in our phenomenologi- 
cal waveform model below which assumes A is the EOS 
dependent parameter. 



VI. PHENOMENOLOGICAL BHNS 
WAVEFORM 



In Paper I and the section above, we assumed that cor- 
relations between A and the other parameters are negli- 
gible. To test this assumption, we must calculate the 
complete Fisher matrix for all parameters, and this re- 
quires us to evaluate partial derivatives with respect to 
all parameters at a single point. As discussed above, 
this is computationally difficult for hybrid waveforms us- 
ing finite differencing. Another approach is to construct 
an analytic BHNS waveform model with free parame- 
ters that are fit to the hybridized numerical waveforms. 
This allows one to interpolate between the available sim- 
ulations and to evaluate derivatives used in the Fisher 
matrix. 



As found in Section|v](Figs. 13 and 14 1 and in Figs. [3]- 
[5] a BHNS waveform is well approximated by a one- 
parameter deformation from a BBH waveform where 
A = [17]. As shown in Figs.[3j|5j throughout the inspi- 
ral, merger, and ringdown, both the amplitude and phase 
of the Fourier transformed waveform monotonically de- 
crease with respect to a BBH waveform as frequency in- 
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waveform as a modification to a BBH waveform 

/ibhns(M/;0) = W(M/;0)r(M/;0)e iA *( M ^, 

(31) 

where the ratio r(Mf; 9) = 

\hBHNs{Mf;9)\/\hBBn(Mf;9)\ is the amplitude correc- 
tion, A$(M/;0) = $bhns(M/;^) - $bbh(M/;0) is a 
phase correction factor, and the 3 physical parameters 
that we will fit our waveforms to are 9 — {t],xbb\, A}. 
We will then fit the quantities r{Mf; 9) and A$(M/; 9) 
to the 134 hybrid waveforms listed in Table [TT] 

For aligned-spin BBH systems, significant work has 
gone into developing analytic waveforms that include the 
complete inspiral, merger, and ringdown, and are cal- 
ibrated to numerical BBH simulations. Unfortunately, 
there is no single BBH waveform model that has been 
calibrated for the ranges of mass ratio Q € [2,5] and 
black hole spin xbh € [-0.5, 0.75] listed in Table [TTJ Fur- 
thermore, there are differences between the various ana- 
lytic BBH waveforms for values of Q and xbh for which 
they have both been calibrated (See the PhenomC and 
EOB waveforms in Figs. [3j-[5]) . Although small, these 
differences are still a non-negligible fraction of the differ- 
ence between a BBH and BHNS waveform. As a result, 
we will have to separately calibrate our BHNS waveform 
model for each BBH model, and our waveform for very 
small values of A, where intrinsic errors in the analytic 
BBH models may dominate over tidal effects, will likely 
not be accurate. The two BBH models we use are the 
frequency-domain PhenomC |40) waveform and the time- 
domain EOB waveform [?7] discussed in Section III 



A. Fit based on PhenomC BBH approximation 

We fit the corrections to the PhenomC waveform, r and 
A$, to our numerical BHNS waveforms. Although the 
PhenomC waveform is not calibrated using BBH wave- 
forms with mass ratios of Q — 5, we will fit r and A<& to 
Q = 5 BHNS waveforms anyway. 



FIG. 14: Same as Fig. |13[ except the hybrid waveforms 
are generated by matching the numerical BHNS waveforms 
to PhenomC inspiral waveforms with the tidal correction ipx- 
For all panels, the matching window has width AM/ matc h = 
0.008 and is centered on M/ mid = 0.022. 



creases and as A increases 3 . We thus write the BHNS 



3 For mass ratios of Q = 4 and 5, the PhenomC amplitude is 
sometimes slightly less than the BHNS waveform with the softest 
EOS during ringdown around Mf ~ 0.1 as seen in Figs.|3]and|5] 
The amplitude of the corresponding EOB waveform, however, is 
always greater than the BHNS amplitudes. 



1. Amplitude fit 

During the inspiral, because parameter estimation is 
much more sensitive to the fractional change in phase of 
the waveform than to the fractional change in amplitude, 
we ignore the very small amplitude correction from tidal 
interactions. During the merger and ringdown, however, 
amplitude corrections are important. We therefore write 
the amplitude correction as 

<Mf,0)= I 1 ab(m f 8) Mf ~ MfA > (32) 

where MJa is the boundary, chosen below, between the 
inspiral and merger for the amplitude fit. We have ex- 
tracted the quantity r/A because, as rj — > (extreme mass 
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ratio limit) or A — > (no matter limit), the waveform 
should approach that of a BBH waveform. We now im- 
pose two requirements on the function B(Mf; 9). (i) The 
amplitude must be continuous at the frequency MfA, so 
B(Mf A ; 0) = 0. (ii) Because the amplitude of the BHNS 
waveform is almost always less than that of the corre- 
sponding BBH waveform, we require B(Mf;8) > for 
Mf > Mf 'a and for all physical values of the parameters: 
rje [0,0.25], xbh € [-1,1], and A > 0. 

Given the above restrictions, we find that a useful fit- 
ting function for the amplitude correction is B(Mf; 9) = 
C(Mf — MfA) D , where C and D are free parameters, 
and MfA — 0.01. With this ansatz, we then do a non- 
linear least-squares fit to determine the parameters C 
and D. We find that over the 134 simulations, D has 
a mean and standard deviation of D ~ 3 ± 0.5, and 
because these parameters are highly correlated, we fix 
D = 3 so that B = C{Mf - Mf A f- We then fit 
each waveform with the single parameter C. The pa- 
rameter C is then fit to the physical parameters. We 
find that for fixed rj and xbh, C is approximately a 
linear function of A. We therefore use the function 
C(?7,xbh,A) = e fco+bl?7+b2XBH + Ae c ° +C1 " +C2 X BH , where 
the parameters {&0; b±, b 2 , Cq, c±, C2} are found with a non- 
linear least-squares fit. We note that this function is pos- 
itive for all physical values of the parameters 77, xbh 5 and 



A. Given the small number of samples for 77 and xbh as 
well as the difficulty in extracting the small tidal contri- 
bution from numerical simulations, we have used as few 
parameters in our fit as possible rather than to over-fit 
noisy data with a large number of parameters. The final 
form of B is 

B{Mf;9) = ( e b °+ b i"+ fe «BH + Ae C0+cir > +C2XBH ) {Mf-Mf A f, 

(33) 

and the best-fit parameters are {60, b\, 62, Cq, c\, c 2 } = 
{-64.985, -2521.8, 555.17, -8.8093, 30.533, 0.64960}. 
Using this fit, we find typical fractional errors in C of 
~ 30% for A > 500. However, errors can be significantly 
larger for A < 500 and Q — 4 and 5, where r is small 
due to the small tidal interaction, and the error is 
dominated by numerical noise and uncertainty in the 
BBH waveform. This is not significant because, as we 
will find, the systematic error that results from poorly 
fitting r is still smaller than the statistical error in A for 
small values of A. 



2. Phase fit 

For the phase of the phenomenological waveform we 
choose the following ansatz 



A$(M/;t 



MMf;9) Mf<MU 
-r,AE(Mf; 9) + M M U; &) + (Mf ~ M/«)^(M/«; 9) Mf > M/« 



(34) 



where "0T is the frequency-domain tidal phase correction 
for the inspiral, and a ' denotes a derivative with respect 
to Mf. In this paper we will use the 1PN accurate Tay- 
lorF2 tidal correction (Eq. Q) for the inspiral. This 
expression explicitly breaks the EOS-dependent contri- 
bution to the phase into three pieces: (i) the contribu- 
tion due to the inspiral tidal correction tp T (Mf;9), (ii) 
the contribution due to the merger-ringdown dynamics 
—r)AE(Mf; 9) which we will fit to numerical simulations, 
and (iii) the term Vt(M/$; §) + (M f - M U)ip' T (M fa; 9) 
that grows linearly after the transition frequency M/$ 
and results from matching the phase and derivative of 
the merger to the tidally corrected inspiral as discussed 
in Section ITOl 

As in the amplitude fit, we have explicitly pulled out 
the quantity r/A in the first term —r/AE because the phase 
of the BHNS waveform should approach that of a BBH 
waveform as 77 — > or A — > 0. We further require the re- 
maining function E(Mf; 9) to satisfy the following con- 
ditions: (i) E(MU;9) = 0, (ii) E'(MU;9) = 0, and (iii) 
E(Mf; 9)>0 for Mf > Mf® and for all physical values 



°f Vi Xbh, and A. In this way, the function E(Mf;9) is 
determined fully by the numerical waveform and is inde- 
pendent of the inspiral tidal term tpT- 

A key feature of this ansatz for the phenomenological 
waveform is that, if we do not change the hybridization 
matching window (M f il M ff), an improved inspiral tidal 
phase term ipT can be swapped in to Eq. (34) without 
requiring one to redo the following fit for E(M f; 9). This 
is useful for estimating how an improved inspiral tidal 
correction effects the measurability of tidal parameters 
for the complete IMR waveform. We note, however, that 
an improved inspiral tidal term will lead to a slightly 
different optimal matching window (Mfi,Mff) for the 
hybridization procedure, and using the optimal matching 
window for the hybridization with the improved inspiral 
tidal term will require one to redo the fit for E. 

We find that each waveform can be accurately fit with 
a function of the form E = G(Mf - Mf 9 ) H , where G 
and H are free parameters, and unlike the amplitude 
fit where Mf A = 0.01, we choose M/$ = 0.02 for the 
transition frequency because it is close to the midpoint 
(M/ m id = 0.022) of the hybrid matching interval. For the 
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134 BHNS waveforms the best fit for the parameter H has 
a relatively narrow range of ~ 2 ± 0.5. This is consistent 
with the leading frequency dependence (ip T oc (M/) 5 / 3 ) 
of the tidal correction in Eq. Q. In addition, the free 
parameters G and H in this fit are highly correlated. We 
thus rewrite E = G(M f — M /o) 5 ' 3 and fit each waveform 
with the single parameter G. For fixed r\ and Xbh 5 G is 
a roughly constant function of A, so we use the following 
form G(?7,xbh) = e 9'>+9^+92XBH+g 3 vxBH _ Xhe f unc ti on 
E can then be written 

E(Mf]6) = e So+Si'7+S2XBH+33TOBH( M J_ M J )5/3 ; ( 35 ) 

where the best-fit parameters are {g$, g\, g2, 53} = 
{-1.9051, 15.564, -0.41109, 5.7044}, and as with the am- 
plitude fit, this parametrization is well defined for all 
possible values of r], xbh, and A. We find that typical 
fractional errors in the fit for G are ~ 30% for A > 500, 
but can be larger for smaller values of A. As with the 
amplitude fit, the large fitting error for A < 500 is not 
significant because the systematic error resulting from 
the poor fit will still be less than the statistical error in 
A. 



B. Fit based on EOB BBH approximation 

Unlike the PhenomC waveform, EOB waveforms are 
not yet available for spins of xbh = 0.75, so we use only 
the 90 waveforms that have — 0.5 < xbh < 0.5 when cal- 
ibrating the fit. As discussed in Section [ill B| we Fourier 
transform the time-domain EOB waveform, then add the 
TaylorF2 tidal correction (Eq. Q) to the phase, and then 
generate a hybrid waveform. We then produce an ana- 
lytic fit to these hybrid waveforms using the same pro- 
cedure and matching parameters (AM/ matc h = 0.008, 
M/mjd = 0.022, Ms, = Mf u and Ms f = M Si + 0.001) 
as with the PhenomC waveforms. 

For the amplitude fit we use for Eq. (33), M/a = 0.01 
and obtain the following coefficients: 

{&0>&1,&2> Ob Cx,C 2 } 

= {-1424.2, 6423.4, 0.84203, -9.7628, 33.939, 1.0971}. 
Typical errors in C are about the same as for the Phe- 
nomC fit. 

For the phase fit we use for Eq. (35), M/$ = 0.02 and 
obtain 

{50, 91, 92, 33} = {-4.6339, 27.719, 10.268, -41.741}. 
Typical errors in G are again about the same as for the 
PhenomC fit. 



VII. MEASURABILITY OF A 

A. Statistical error 

Using the analytic BHNS waveform based on the Phe- 
nomC BBH waveform developed in the previous sec- 
tion, we can now evaluate the Fisher matrix for a sin- 
gle gravitational-wave detector using the complete set of 



waveform parameters {lnD c g, f\t c , 4> c , hiM, In 77, xbh, 
A 1 / 5 }, where f\ is some fiducial frequency such as 1 Hz, 
and as in Paper I we use A 1 / 5 because it is approximately 
proportional to the more familiar NS radius. We have cal- 
culated the 1-a uncertainty in A 1 / 5 for both the broad- 
band aLIGO [57] and ET-D detector configurations [58] 
shown in Fig. |11| Errors are shown in Figs. [15] and [16 for 
broadband aLIGO and ET-D respectively, and are scaled 
to an effective distance of 100 Mpc as was done in Pa- 
per I. We note that the results here for the Q = 2 and 3 
nonspinning waveforms are similar to those presented in 
Figs. 12 and 13 of Paper I. This indicates that coherently 
adding the inspiral tidal interactions to the merger and 
ringdown, and considering correlations between A and 
the other parameters roughly cancel each other. 

There are several trends to notice in the uncertainty 
. In general, <7 A i/s increases with increasing mass 
ratio Q. This is not surprising since the inspiral tidal con- 
tribution tpx (Eq. (R) to the waveform phase, which has 
a significant impact on the measurability of A both before 
and after the inspiral-merger transition, is a decreasing 
function of the mass ratio. In addition, the amount of 
tidal disruption before the plunge, as well as its imprint 
on the waveform, decreases as the mass ratio increases. 
However, there are two competing effects that help to 
minimize the increase in uncertainty er A i/5 as Q increases. 
First, the amplitude during the inspiral which to Newto- 
nian order scales as \h(f)\ oc A4 5 ^ 6 f~ 7 ^ 6 /D e g increases 
as the mass ratio increases for a fixed NS mass. Sec- 
ond, for higher mass ratios, the EOS dependent merger 
dynamics occur at lower frequencies, closer to the mini- 
mum of the noise PSD (Fig. 11 ). 

On the other hand, cr A i/s decreases with increasing 
black hole spin xbh ■ This effect can be understood from 
Fig. [7] where the amplitude and phase difference between 
a BHNS waveform and BBH waveform with the same 
parameters increases as the BH spin increases, and the 
amplitude cutoff occurs at a lower frequency where the 
detector is more sensitive. Physically, we expect the EOS 
dependence to be greater for higher spins because the 
BH ISCO decreases with spin, allowing the NS to be- 
come more tidally distorted before passing through the 
ISCO and plunging into the black hole. In addition, for 
high spins, the orbital decay will take longer because 
the system must radiate away sufficient angular momen- 
tum for the final Kerr black hole to have spin parameter 
Xbh' < !• As a result, the waveform has more chance 
to deviate from a BBH waveform towards the end of the 
inspiral and into merger. 

We also note that in general, cx A i/5 decreases as a 
function of A 1 / 5 . This occurs because the departure 
from BBH behavior, given by Eqs. p2| and (34), is a 
strongly increasing functions of A 1 / 5 . However, as seen in 
Fig.[l5]for broadband aLIGO, we find that for Q = 2 and 
Xbh ^ 0.5, the error <r A i/5 begins to increase again for 
large A 1//5 . Because the Fisher matrix element r A i/s A i/5 
is a monotonically increasing function of A 1 / 5 , the in- 
crease in <7 A i/5 is, therefore, due to an increase in the 
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covariance with the other parameters for this particu- 
lar set of parameters and noise curve. Finally, we find 
that although the uncertainty <t A i/s for the parameter 
values (Q = 2, X bh = 0, M NS = 1.35M©), (Q = 3, 
Xbh = 0.5, M NS = 1.35M ), and (Q = 5, xbh = 0.75, 
M NS = 1.35M ) only varied by ~ 50% in Fig. [IIJ they 
vary by ~ 100% in Figs. [15] and [16] This again results 
from accounting for the covariance with the other param- 
eters. 



B. Systematic error 

The analytic BHNS waveform developed in Section |Vl| 
does not exactly match the hybrid waveforms that it is 
calibrated against. This leads to a bias in estimating 
the BHNS parameters. In this subsection, we estimate 
the systematic error that results from using our analytic 
waveform instead of the hybrid waveforms as templates, 
then we use it as a criteria to determine the quality of 
the analytic fit. 

In section |IV[ we gave an approximate expression for 
the systematic error 



A8. 



syst 



-T. 1 (5h(9 A ),d 3 h A A )), (36) 



where Sh(t;6) — tiA(t;9) — h,E(t;9) is the difference be- 
tween the approximate fit and exact waveform. Here, we 
refer to the hybrid BHNS waveform as the exact wave- 
form /ig(t;0), the analytic fit as the approximate wave- 
form h,A(t;6), and 6 a is the best estimate of the true 
parameters using the approximate fit as the template. 
The waveform derivatives and Fisher matrix are then 
calculated from the analytic waveform fit as was done for 
the statistical error. 

An analytic waveform is useful for parameter estima- 
tion when the systematic error is a small fraction of the 
statistical error for a given gravitational wave detector. 
We will use as our criteria for a sufficiently accurate fit 
to the numerical waveform the requirement that the sys- 
tematic error in each parameter be less than the statis- 
tical error for an optimally oriented BHNS system ob- 
served at 100 Mpc. For larger distances or a less sensi- 
tive detector, the systematic error will be a smaller frac- 
tion of the total error. For a waveform with the form 
h(t) — jy^g(t) and a PSD with overall amplitude fac- 
tor A n such that S n (f) = A^ l R n (f), the statistical errors 
<7j scale as <7j oc A n D c tf. On the other hand, as can be 
found from Eq. p6| , the systematic error is independent 
of both A n and £> c ff • The ratio of systematic to statistical 
error is therefore 



A0: 



syst 



1 



A n D. 



off 



(37) 



Because the PSD for broadband aLIGO and ET-D have 
roughly the same shape, but differ in amplitude A n by 
a factor of ~ 10, we expect the systematic to statistical 



error ratio to differ by about a factor of 10, and therefore 
the systematic error will be far more important for ET. 

We find, for the BHNS fit based on the PhcnomC 
waveform, the ratio of systematic error to statistical 
error in A 1 / 5 for broadband aLIGO for a binary at 
D cB = 100 Mpc is |AA^ s 5 t |/o- A i/5 ~ 0.15 ±0.15 with 
a maximum value of 0.65 for the 134 waveforms listed 
in Table [H] However, the systematic error is strongly bi- 
ased by the Q = 5 waveforms for which the PhenomC 
waveform has not been calibrated to BBH simulations 
(Table [Tl]). For the 101 simulations with Q < 4, we find 
|AA^|/tr A i/ 8 ~ 0.09 ±0.08 with a maximum value of 
0.51, and we also note that this ratio is > 0.2 only when 
Q>4. 

For ET-D, the systematic error is roughly the same as 
for broadband aLIGO, and the ratio of systematic to sta- 
tistical error is about an order of magnitude larger as ex- 
pected. Specifically, at D cS — 100 Mpc, |AAgyg t |/cr A i/5 ~ 
1.3 ± 1.2 with a maximum value of 6.3 for the 134 wave- 
forms. Only including the 101 waveforms where Q < 4, 
|AA^ s 5 t |/o- A i /5 ~ 1.1 ± 0.9 with a maximum value of 5.3, 
and this ratio is > 2.5 only when Q > 4. We therefore 
conclude that the BHNS fit is sufficient for aLIGO. How- 
ever when the effective distance D c g is less than a few 
hundred Mpc, the systematic error will become compa- 
rable to the statistical error in some cases for ET-D. 

We also calculate the systematic error for the fit based 
on the EOB waveform. Unfortunately, because the EOB 
waveform is time-domain and the waveform enters the 
detector band at around 10 Hz for aLIGO and 1 Hz for 
ET, very long waveforms are needed to evaluate the sta- 
tistical and systematic error. This is possible, but time 
consuming, so we use a more crude estimate of the errors. 
Eqs. (24) and (36) can be approximated by ignoring the 



covariance between A and the other parameters: 



(7 A l/5 



(d A i /5 h A (0A), 0aV5/u(0^))Im/=o°o 



-1/2 



(38) 



(AA 1/,5 ) syst « - 



(Sh(9 A ), d A i /5 h A (0A))\Mf J^oi 



(<9 A i/5 h A (9a), d A i/5 h A {6 A )) Im /= ^° i 



(39) 



iM/=oo 

Im/=o.oi 



IS 



where the integral in the inner product (• 
evaluated from a frequency of Mf — 0.01 to oo. We 
also note that the hybrid waveform Ke and analytic fit 
Ka are identical below Mf = 0.01, so Sh = below this 
frequency. 

We have evaluated these quantities for the BHNS fits 
based on the PhenomC and EOB waveforms. The sta- 
tistical errors are the same to approximately ±10%. The 
systematic error for the EOB fit, however, is typically ~ 2 
times larger than for the PhenomC fit. This is somewhat 
surprising given that the EOB waveform is a much bet- 
ter approximation to BHNS waveforms with the smallest 
values of A. However, we note that most of the adjustable 
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FIG. 15: l-<7 error <7 A i/5 for various values of the mass ratio, BH spin, and tidal deformability. NS mass is fixed at 1.35 




FIG. 16: Same as Fig. |15| but with the ET-D noise curve. Uncertainty o" A i/s is an order of magnitude smaller, 
parameters, such as the matching window, as well as the 



form of the functions r and A$ (Eqs. (32)-(35|) were 
optimized for hybrids based on the PhenomC waveform. 
We therefore do not claim that the PhenomC waveform 
in general produces a better BHNS analytic fit. 



VIII. DISCUSSION 

A. Summary 

We have examined the ability of gravitational wave de- 
tectors to extract information about the EOS from obser- 



21 



vations of BHNS coalescence for black holes with aligned 
spin. In Paper I, we found that the EOS parameter that 
is best measured during the merger and ringdown, for 
systems with nonspinning black holes, is consistent with 
the tidal dcformability A. We have now found that this is 
also true for systems with aligned black hole spins. Fur- 
thermore, coherently joining the tidally corrected inspi- 
ral, which analytically depends on A, to the merger and 
ringdown dramatically improves the alignment of the er- 
ror ellipses with A in Fig. [14] as well as the measurability 
of A by up to a factor of ~ 3 in some cases over just the 
merger and ringdown. 

In order to examine the correlations between A and 
the other parameters, we constructed an analytic IMR 
waveform based on the frequency-domain, aligned-spin 
PhenomC BBH waveform model [ID] as well as the time- 
domain EOB waveform model [47] , and we calibrated this 
waveform model to our hybridized numerical waveforms. 
Although A does correlate with the other parameters, 
the correlations are not nearly as strong as correlations 
between the other parameters. Overall, the correlations 
reduce the measurability of A by approximately a fac- 
tor of 3. The above two effects roughly cancel out and 
the results for <t a i/5 are therefore similar to the results 
presented in Paper I which neglected these two effects. 

In addition, we examined the agreement between the 
hybrid BHNS waveforms and our analytic waveforms. 
The agreement is good enough that systematic errors will 
be smaller than statistical errors for aLIGO. However, for 
ET, systematic errors will matter if BHNS systems are 
observed with effective distances of less than a few hun- 
dred Mpc. 

B. Future Work 

There is currently enough uncertainty in the modeling 
of BBH systems and inspiral tidal interactions, as well as 
in BHNS simulations and the hybridization procedure, 
that the analytic BHNS waveforms presented here can 
only be considered preliminary. Additional systematic 
errors exist. The waveform model has the ability to in- 
corporate improvements in the inspiral point-particle and 
tidal interactions. However, incorporating EOS depen- 
dent corrections into the merger-ringdown waveform re- 
quires an accurate understanding of the late dynamics for 



both BBH and BHNS systems. This analytic waveform 
fit will therefore need to be re-calibrated when improve- 
ments are made to the late dynamics. On the numer- 
ical side for BHNS systems, key improvements beyond 
the standard issue of convergence, would be to minimize 
eccentricity in the numerical waveforms and to increase 
the number of orbits before merger so that the numer- 
ical waveform can be matched to the inspiral waveform 
at lower frequencies where errors in the analytic point- 
particle and tidal interactions are smaller. On the an- 
alytic side, we are restricted to analyzing correlations 
between parameters for systems with small mass ratios 
and moderate black hole spins because IMR BBH mod- 
els have not yet been calibrated to BBH simulations with 
larger mass ratios and spins. 

The analytic waveform models presented here, and in 
particular the one based on the frequency-domain Phe- 
nomC waveform, can be, without too much difficulty, in- 
corporated into Markov Chain Monte Carlo and Nested 
Sampling algorithms used for Bayesian parameter esti- 
mation for networks of gravitational-wave detectors. A 
full Bayesian analysis will then make it possible to assess 
the true nature of the statistical and systematic errors 
beyond the Fisher matrix approximation, by injecting hy- 
brid BHNS waveforms into detector noise and attempting 
to recover their parameters with the analytic waveform 
template. 
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